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ABSTRACT 

The gravitational lensing properties of cosmological halos depend upon the mass distri- 
bution within each halo. The description of halos as nonsingular, truncated isothermal 
spheres, a particular solution of the isothermal Lane-Emden equation (suitably mod- 
ified for A / 0), has proved to be a useful approximation for the halos which form 
from realistic initial conditions in a CDM universe. We derive here the basic lensing 
properties of such halos, including the image separation, magnification, shear, and 
time-delay. We also provide analytical expressions for the critical curves and caustics. 
We show how the scale-free results we derive yield scale-dependent lensing properties 
which depend upon the cosmological background universe and the mass and collapse 
redshift of the lensing halos, according to the truncated isothermal sphere (TIS) model 
of CDM halos derived elsewhere. We briefly describe the application of these results 
to the currently-favored ACDM universe. 

Key words: cosmology: theory - dark matter - galaxies: clusters: general - galaxies: 
formation - galaxies: halos - gravitational lensing 



1 INTRODUCTION 



Turner, Ostriker, 



The grplLtiLloxial leiiHhig of dlHlauL aonreea haa in recent 
years beiuine unn uf ftfac musL ijuw aiful Luuls in uuMiiva- 

and lef- 



^ ■ tional cosmology (see, foi example, Lioiicail 2001 



Gott 1984; Hinshaw & Krauss 1987 



1995 ; Premadi et al. 2001 



Narayan fc Whi te 1988]: " Blandford et al. 1991 ; Jaroszynski 



1991] , |1992|; |Kochanek 1995t |Premadi, Martcl. fc Matzner 



erences therein). Since the effects of gravitational lensing 
depend upon the redshift of the source, the cosmological 
background, and the distribution of matter in the universe, 
they can be used to constrain the cosmological parameters 
and the primordial power spectrum of density fluctuations 
from which structure originates. In addition, many of the 
effects produced by gravitational lenses, such as image mul- 
tiplicity, separations, and time delay, depend strongly upon 
the maUei disLiibuliuu inside Lire lenses. Hence, measuie- 



HRusin fc Ma 200 ip . As the theory 
of halo formation in the CDM model has advanced in recent 
years, however, the halo mass profiles adopted for lensing 
models have been refined to reflect this theory. Numerical 
simulations of large-scale structure formation in Cold Dark 
Matter (CDM) universes predict that galaxies and clusters 
have a singular density profile which approaches a power law 
p oc r~ n at the center, with the exponent n ranging from 1 



ments c f Lhcjc effect j can piuuidc a iniiuuc luul fui piubing 
the ma tH-^fcsfcrilB'fffcjemnsiete'eeHapBeeHel^ 

and clu sters, providing the only direct measurement of their 



1997|;[Tbr- 
1997, pO 



to 1.5 (ICole fc Lacey 1996|; [Navarro, Frenk, fc White 199( 



men, Bouchet, fc White 1997|; jFukushige fc Makino 
001a|jb|; [Moore et al. 19981 |1999|; [Huss, Jain, fc Stein 



al. 2000 



metz 1999E |Ghigna et al. 200(!| ; [ling fc Suto 2000} |Klypin~ 



Power et al. 2002). These results are in apparent 



dark matter content, and constraining the theory of their 
formation and evolution. 

Until recently, the internal structure of halos adopted 
in lensing studies was generally some gravitational equi- 
librium distribution, either singular or nonsingular (e.g., 
King model, singular isothermal sphere, pseudo-isothermal 
sphere), not necessarily motivated direc tly by the theory of 
cosmological halo fomialioii (s^e, e.g., Young el al. 1980[ ~ 



conflict with observations of rotation curves of dark-matter 
dominated dwarf galaxies and low surf ace brightness galax- 
ies, which favor a flat-density core (cf . [Primack et al. 199!: 



Burkert fc Silk 1999; Moore et al. 1999; |Moore 2001). On the 



scale of clusters of galaxies, observations of strong gravita- 
tional lensing of background galaxies by foreground clusters 
also favor the presence of a finite-density core in the cen- 



ters o f clusters (see, e.g., Tyson, Kochanski, 
1998;). 
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Several possible explanations have been suggested in 
order to explain this discrepancy. The rotation curve data 
might lack sufficient spatial resolution near the center to 
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distinguish unambiguously between a density profile wi th a 
flat-de nsity core and one with a singular profile (e.g. van 
den Bo:ich & Swaters 2001). Attempts have also been made 
to improve the numerical resolving power of the simulations 
to obtain a more accurate determination of th e slope of the 
pre dicted density profil es at small radii (e.g 
199ST 



Moore et al 



Power et al. 2002). However, if the flat-core interpre- 



tation of the observations and the singular cusps predicted 
by the numerical simulations are both correct, then the sim- 
ulation algorithms may be ignoring some physical process 



In particular, the TIS mass profile agrees well with the fit 
by NFW to N-body simulations (i.e. fractional deviation 
of ~ 20% or less) at all radii outside of a few TIS core 
radii (i.e. outside a King radius or so). It also predicts the 
internal structure of X-ray clusters found by N-body and 
gasdynamical simulations of cluster formation in the CDM 
model. For example, the TIS model reproduces to great ac- 
curacy the mass-temperature and radius-temperature virial 
relations and integrated mass profiles de rived empirically 



from the simulat ions of cluster formation (Evrard, Metzler 



which would, if included, serve to flatten the halo density & Navarro 1996). The TIS model also successfully repro- 



profiles at small radii relative to the results for purely grav- 
itational, N-body dynamics of cold, collisionless dark mat- 
ter, while retaining the more successful aspects of the C DM 
model. For example, gasdynamic al processes (see, e.g. El- 
Zant, S llosman, & Hoffman 2001) and a modification of the 



microscopic properties of CDM, such as the proposal of self- 
interacting dark matter (Spergel & Steinhardt 2000), both 



have the potential to lower the central density of halos and 
possibly reconcile simulations with observations. 

Lensing by the two kinds of halo mass profiles, singular 
versus flat-core, will be different. This has led to attempts to 



predict the differences expected if the halos have the singu- 
lar cusp of th e NFW or Moore profiles or else a profile with a 
flat cor e (e. g. Kochanck 1995 ; Keeton fc Madau 2001 ; Rusin 
& Ma iool ; jWyithe, Turner, fc Spcrgcl 200 1| ; |Takahashi & 



duces to high precision the mass-velocity dispersion relation 
for clusters in CDM si mulations of the H ubble volume by 
the Virgo Consortium (Evrard et al. 2002), including its de- 
pendence on redshift for different background cosmologies. 
The TIS model also correctly predicts the average value of 
the virial ratio in N-body simulations of halo formation in 
CDM. 

The TIS profile matches the observed mass profiles of 
dark-matter-dominated dwarf galaxies. The observed rota- 
tion curves of dwarf galaxies are generally well fit b y a den- 
sity p rofile with a finite density core suggested by Burkert 
(1995D , given by 



p{r) = 



Po,s 



(r/r c + l)(r 2 /^ + l) 



(1) 



Chiba ^001) |Li fc Ostriker 2002; ) . Singular profiles like that of 
NFW are physically motivated by the N-body simulations, 
and the latter have been used to place these halo profiles 
empirically in a proper cosmological context which permits 
statistical predictions for the CDM model. The nonsingu- 
lar profiles which have been adopted to contrast with these 
singular ones, however, are generally no more than param- 
eterized, mathematical fitting formulae, with no particular 
physical model to motivate them or put them in a proper 
cosmological context. 

We have developed an analytical model for the post- 
collapse equilibrium structure of virialized objects that 
condense out of a cosmological background universe, ei- 
t her matter-dominated or flat with a cosmological c onstant 
( Shapiro, Iliev, fc Raga 199!: , hereafter Paper I; Qiev & 
Shapirc 2001a , hereafter Paper II). This Truncated Isother- 
mal Sphere, or TIS, model assumes that cosmological halos 



The TIS model gives a nearly perfect fit to this profile, with 
best fit parameters Po,b/po,tis = 1-216, r c /ro,Tis = 3.134, 
correctly predicting the maximum rotation velocity f max 
and the radius r max at which it occurs. The TIS model 
can also expl ain the mass profile with a flat density core 
measured by Tyson, Kochanski, & DelP Antonio (1998) for 



cluster CL 0024+1654 at z = 0.39, using the strong grav- 
itational lensing of background gal axies by the cluster t o 
infer the cluster mass distribution ( fchapiro fc Iliev 200 1[ ). 
The TIS model not only provides a good fit to the projected 
surface mass density distribution of this cluster within the 
arcs, but also predicts the overall mass, and a cluster veloc- 
ity dispersion in cl ose agreement with t he value a v = 1150 
km/s measured by Dressier et al. (1999). 



Several authors have studied the effect of lens i ng by 



halos with a f lat-density core ( Jaroszyhski 1991, 199S 
Kochanek 1995 ; Premadi, Martel, & Matzner 1998 ; Premadi 



form fr om the collapse and virialization of "top-hat" density et al. 2001| ) or by NFW or Moore profiles that have been gen 



perturbations and are spherical, isotropic, and isothermal. 
This leads to a unique, nonsingular TIS, a particular solu- 



tion of the Lane-Emden equation (suitably modified when 
A 7^ 0). The size r t and velocity dispersion ov are unique 
functions of the mass M and formation redshift z co l\ of the 
object for a given background universe. The TIS density pro- 
file flattens to a constant central value, po, which is roughly 
proportional to the critical density of the universe at the 
epoch of collapse, with a small core radius ro ~ n/30 (where 
o~y — 4irGporQ and rp _ ruing /3, for th e "King radius" 



eralized, so that the inner slope of the density pro f ile is ar- 
bitrary ([Keeton fc Madau 200 1|; |Rusin fc Ma 2001] ; [Wyithe 



Turner, fc Spergel 2001| ; [Li fc Ostriker 2002j ). These par- 



king, defined by Binney & Tremaine 1987, p. 228) 



Even though the TIS model does not produce the cen- 
tral cusp in the density profile of halos predicted by nu- 
merical CDM simulations at very small radii, it does repro- 
duces many of the average properties of these halos quite 
well, suggesting that it is a useful approximation for the 
halos whic h result from more realis tic initial conditions (Pa- 
pers I, II; Iliev & Shapiro 2001b and references therein). 



ticular density profiles are essentially mathematical conve- 
niences without physical motivation. There is no underlying 
theoretical model in these cases that was used to predict the 
value of the core radius or the departure of the inner slope 
of the density profile from the value found by N-body sim- 
ulations of CDM. By contrast, the TIS model is based on 
a set of physical assumptions concerning the origin, evolu- 
tion, and equilibrium structure of halos in CDM universes. 
Observations of gravitational lenses have the potential to 
distinguish between the TIS profile and singular ones like 
the NFW profile, as several observable properties of gravi- 
tational lenses will be strongly affected by the presence, or 
absence of a central cusp in the density profile. One exam- 
ple of an important observable that can distinguish between 
various density profiles is the parity of the number of im- 
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ages. Lenses with nonsingular density profiles, such as the 
TIS, obey the odd number theorem. The number of images 
of a given source is always odd, unless the source is extended 



Falco 1992 



and saddles a caustic (see Schneider, Ehlers, 
hereafter SEF, p. 172). Lenses with singular profiles, like the 
singular isothermal sphere, the NFW profile, or the Moore 
profile, need not obey this theorem, even for point sources. 
Most observed multiple-image gravitational lenses have ei- 
ther 2 or 4 im ages, and this may argue against profiles with 
a central core (Elusin & Ma 2001). There are, however, other 
possible explanations for the absence of a third or fifth im- 
age. That image tends to be very close to the optical axis, 
and might be hidden behind the lens itself. Also, it is usually 
highly demagnified, and might be too faint to be seen. 

We can use the TIS solution to model observed gravita- 
tional lenses individually. Alternatively, we can use the ob- 
servations collectively to constrain the distribution of halo 
properties as characterized by the TIS solution. These prop- 
erties, core radius, velocity dispersion, central density, and 
so on, depend upon the mass of the lensing halos and the 
redshift at which they form. Observational constraints on 
the statistical distribution of these properties will, in turn, 
impose constrains on the cosmological parameters and the 
primordial power spectrum of density fluctuations. 

The problem of studying gravitational lensing of distant 
sources in an inhomogeneous universe can be divided into 
two parts. The first part consists of determining the intrinsic 
properties of the lenses. In particular, we need to determine 
the relationship between the observables (image multiplic- 
ity, magnification, brightness ratio, sheer, image separation, 
time delay, . . .) and the lens parameters. The second part 
consists of tying the lens properties to the cosmology. This 
involves using cosmological models of structure formation 
and evolution to determine the statistical distribution of the 
lens parameters, the clustering properties of the lenses, and 
the nature of their environments. In this paper, we focus on 
the first part, determining the intrinsic properties of individ- 
ual lenses, which is an essential building block. The second 
part will be the subject of forthcoming papers. 

The remainder of this paper is organized as follows. In 
§2, we derive the lens equation. In §3, we compute the critical 
curves and caustics. In §4, we study the properties of multi- 
ple images: separation, magnification, brightness ratios, and 
time delay. In §5, we place the scale-free description of these 
properties in the cosmological context of the CDM model 
and explain how the dimensionless parameters of §2— §4 are 
related by the TIS model to the properties of lensing halos 
in physical units. Summary and conclusion are presented in 
§6. 



2 THE LENS EQUATION 

Figure 1 illustrates the lensing geometry. r\ and £ are the 
position of the source on the source plane and the image 
on the image plane, respectively, a is the deflection angle, 
and Dl, Ds, and Dls are the angular diameter distances 
between observer and lens, observer and source, and lens 
and source, respectively. The lens equation is 
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(2) 
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Figure 1. The lensing geometry: the dots indicate the location 
of the observer, lensing galaxy, and source. r c is the core radius 
of the galaxy, and 77 is the distance between the source and the 
optical axis. The angular diameter distances D^, D^g, and Dg 
are also indicated. 



Notice that since the lens is axially symmetric, we can write 
the quantities 77, £, and a as scalars instead of 2-component 
vectors. We introduce a characteristic length scale £0, and 
rescale the positions and deflection angle, as follows: 



y 



D L rj 



Dsta 

.-f, 

DlDlsol 
ot = — — • 

The lens equation reduces to 
y = x — a(x) . 



(3) 
(4) 
(5) 

(6) 



To compute the deflection angle a, we first need an expres- 
sion for the projected surface density ct(£). The density pro- 
file of the TIS is well-fitted by the following approximation: 



p{r) 



Po 



A 



B 



+ r 2 /rg b 2 +r 2 /rf ) 



(7) 



(Paper I) where A = 21.38, B = 19.81, a = 3.01, b = 3.82. 
The projected surface density is given by 



cr(0 



p{r)dz, 



(8) 



where £ is the projected distance from the center of the TIS, 
and 2 = (r 2 — £ 2 ) x / 2 _ This expression assumes that the fit 
given by equation (Q) is valid all the way to r = 00. Actually, 
it is only valid up to a truncation radius r t ~ 30ro (the ac- 
tual value is 29.4 for an Einstein-de Sitter universe [Paper I], 
slightly different for an open matter-dominated universe, or 
a flat universe with a cosmological constant [Paper II] ) . One 
could always set the limits in equation (^) to ±(r 2 — £ 2 ) 1//2 . 
However, it turns out that the change in a would be small, 
simply because the value of 30 is significantly larger than 
both a and b. Also, the density p(r) is nonzero at r > 30ro, 
and ignoring its contribution entirely would be incorrect. 
For the sake of simplicity, we shall assume that equation (^) 
remains a good approximation out to r = 00. We substitute 
equation (M) in equation (H), and get 



A 



B 



{a?rl+e) 1/2 (& 2 rg + e 2 ) 1/2 



(9) 
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[This result was also d erived by Nfatarajan fc Lynden-Bell 
(1997^Jind Iliev (200C)]. For spherically symmetric lenses, 
the deflection angle is given by 



a(x) 



x'^-dx' 

Ocrit 



(10) 



[SEF, eq. (8.3)], where cr cr it is the critical surface density, 
given by 

c 2 Ds f . u 
acrit = 4,GD L D LS ' 

where c and G are the speed of light and the gravitational 
constant, respectively, and Dl, Ds, and Dls are the angu- 
lar diameter distances between observer and lens, observer 
and source, and lens and source, respectively. We substi- 
tute equation ([]) into equation (|lo|), eliminate £ using equa- 
tion Ql), and set the characteristic scale £o equal to ro- We 
get 

= 2wo^ 2+x2 1/2 
Merita; L 

-B{b 2 +x 2 ) 1/2 - Aa + Bb] . (12) 

We now introduce the dimensionless central surface density, 
or central convergence, k c , defined by 

Ccrit CTcrit \ « J 

and use this expression to eliminate <7 cr it in equation (|l2|). 
It reduces to 



2abn c 



(Ab - Ba) 



r a i 2 . 2n 1/2 

■ [A(a + x ) ' 



-B{b 2 +x 2 ) 1/2 - Aa + Bb] 
This expression has the following limiting cases: 
k c x , x <C a, b; 

x ^> a, b . 



a(x) 



2ab(A - B)k c 



(14) 



(15) 



L Ab- Ba 
The final form of the lens equation is 
2abn. 



V = x- 



(Ab- Ba)x 



2sl/2 



[A(a 2 + x 
-B{b 2 +x 2 ) 1/2 -Aa + Bb] 



(16) 



[This result was also obtained by Chiba & Takahashi (2001) 



3 CRITICAL CURVES AND CAUSTICS 
3.1 Solutions 

The determination of the critical curves is quite trivial for 
axially symmetric lenses. The dimensionless interior mass 
m(x) is related to the deflection angle a(x) by 

i \ — i \ 2abn c r 2 . 2ni/2 

m(x) = a(x)x = Ba j [ A ( a + x ) 

- B{b 2 + x 2 ) 1/2 - Aa + Bb] . (17) 
[SEF, eq. (8.3)]. Tangential critical curves are defined by 
m(x) 2abn c 



/2 



(Ab - Ba)x- 



• [A(a 2 + x 2 ) 1 
- B(b 2 + x 2 ) 1/2 -Aa + Bb] =1. (18) 
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Figure 2. Top: m(x)/x 2 versus x, for 3 particular values of n c . 
Tangential critical curves, denned by m(x)/x 2 = 1, can only occur 
for k c > 1. Bottom: d(m/x)/dx versus x, for 3 particular values 
of fc c . Radial critical curves, defined by d(m/x)/dx = 1, can only 
occur for k c > 1. 



This can be turned into a fourth-degree equation for x 2 . 
Even though such equation can be solved analytically, this 
is a case where a numerical solution is preferable. But first 
let us investigate the existence of a solution. According to 
equations f jl5p and (0), m(x) = k c x 2 at small x. As x in- 
creases, the dependence of m(x) on x drops from m(x) oc x 2 
to m(x) oc x. If k c < 1, m(x) "starts up" below x 2 and 
can never raise above it, and therefore equation (jl8|) has no 
solution. If k c = 1, equation (jl8|) has a unique solution at 
x = 0. Finally, if k c > 1, there is a unique solution x > 0. 
This is illustrated in the top panel of Figure 2. This result 
was expected: critical curves are associated with the phe- 
nomenon of multiple imaging. For axially symmetric lenses 
with monotonically decreasing surface density (da/d£ < 0), 
multiple images can occur only if k c > 1 (SEF, p. 238). 
Radial critical curves are defined by 



d(m/x) 
dx 



2abn r 



(Ab - Ba)x 2 



Aa 2 



+ 



Bb 2 



(ft2 + 3,2)1/2 



(a 2 +X 2 ) 1 / 2 



Aa - Bb 



= 1 . 



(19) 



According to equation (|l5|), d(m/x)/dx = da/dx — k c 
at small x. Equation ( |ig| ) clearly shows that d(m/x)/dx is a 
monotonically decreasing function of x. Hence, equation ( |l9| ) 
has a nonzero solution only if n c > 1. This is illustrated in 
the bottom panel of Figure 2. 

We have solved equations (|l8|) and ( Jis| ) numerically for 
the tangential critical radius xt and radial critical radius x r - 
The solutions are plotted in Figure 3, as functions of k c . Also 
plotted is the radial caustic radius y r , obtained by substi- 
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Figure 3. Radii of the radial critical circle, x r , tangential critical 
circle, xt, and radial caustic, y r , versus k c . 



tuting the value of x r into equation (|l6|). (The value of y r 
we obtain is actually negative, but this reflects the fact that 
the equation locates the point in the plane of the sky which 
is the intersection between the radial caustic circle projected 
onto that plane and a line in this plane from some point on 
the radial critical circle through the origin at the center of 
the lens; the intersection point on the radial caustic circle is 
on the opposite side of the origin. The actual radius of the 
caustic circle, then, is the absolute value of y r in that case.) 
Both xt and y T increase rapidly with k c , while the value of 
x r levels off. We can easily find the asymptotic behavior in 
the limit of large k c . To find Xt, we take the limit k c — » oo in 
equation ([l8|). To satisfy this equation, either the factor x 2 
in the denominator must diverge, or the quantity in bracket 
must vanish. But this quantity is zero only for x — 0, and 
approaches as x 2 , cancelling the x 2 in the denominator. 
Hence, the only way to satisfy equation ( p^ ) in the limit 
tic — > oo is by having x — > oo in the denominator. We take 
the limit x — > oo in equation (|l^), and keep the first and 
second leading terms. We get 



2abn c 
Ab- Ba 



B + 



Bb - Aa 



(20) 



Setting x = oo in the last term, we get x « 2ab(A — 
B)n c /(Ab — Ba). We then substitute this value in the last 
term of equation (pp|), and get 



x t 



2ab(A - 



Ab- 



B) Bb 
Ba A ■ 

= 1.638« c 



Aa 
~B~ 
-7.210, 



k c > 1 ■ 



(21) 



To find x r , we take the limit k c — » oo in equation (|l9|). 
In this case, the quantity in brackets first increases with x, 
then decreases and drops to at a finite value x — 6.479. 
Hence, unlike equation ([ll]), it is the bracket, not the factor 
of x 2 in the denominator, that prevents the expression from 
diverging as k c ^ oo. The asymptotic limit is therefore 



Finally, we substitute this value of x r in equation (|l6|), and 
take the limit k — > oo. We get 



y r = 2.425k c - 6.479 , k c > 1 . 



3.2 Illustrative example 



(23) 



x r = 6.479 , 



Kc > 1. 



(22) 



Using a simple ray-tracing algorithm, we computed the im- 
age(s) of a circular source of diameter Ay — 1, created by 
a TIS with central convergence k c = 4.015. The results are 
shown in Figure 4 for 8 different locations of the source, 
ranging from y — 8.0 to y — 0.0. For each case, the left panel 
shows the source and the caustic circle (y r = 5.640) on the 
source plane, and the right panel shows the images(s), the 
radial critical circle (x r = 3.334), and the tangential critical 
circle (xt = 9.783) on the image plane. For the cases y — 8.0 
and 6.0, only one image appears. At y — 5.4, the source 
overlaps the caustic, and a second, radially-oriented image 
appears on the radial critical circle. At y = 4.8, the source 
is entirely inside the caustic, and the second image splits in 
two images, located on opposite sides of the radial critical 
circle, forming with the original image a system of 3 aligned 
images. As the source moves toward y — 0, the central image 
moves toward x = and becomes significantly fainter, while 
the other images move toward the tangential critical circle 
and become bright, elongated arcs. At y = 0, the two arcs 
have merged to form an Einstein ring located on top of the 
tangential critical circle, while the central image, very faint, 
is still visible in the center. 

The central image is always located inside the ra- 
dial critical curve, which has a maximum radius of 6.479 
(eq. |p2[). The density profile of the TIS extends to a cutoff 
radius of order 29.4 (Papers I and II) [Note: We are neglect- 
ing here the effect of matter distributed outside the virial 
radius of our lensing halo.]. Hence the central image is al- 
ways located "inside" the TIS and will be visible only if the 
lens is transparent. Incidently, all the images shown in Fig- 
ure 4 are located inside the TIS, because the value we used 
for k c is rather modest. Setting x t = 29 .4 in equation (pl[), 
we get k c = 17.9. This is the value of k c for which the tan- 
gential critical circle is located along the edge of the TIS. 
For opaque lenses, the image located between the tangential 
and radial critical circles (leftmost images in Fig. 4) might 
be visible only if k c > 17.9. For n c < 17.9, only the out- 
ermost image, located outside the tangential critical circle, 
might be visible. 



4 PROPERTIES OF MULTIPLE IMAGES 
4.1 Image separation 

The locations of the images are computed by solving the lens 
equation (|l6|). After rewriting this equation as x — y = a(x), 
we can solve it graphically. In Figure 5, we show the multiple 
image diagram for the particular case k c = 5.005. The solid 
curves shows a(x), while the dotted lines show x — y for 
particular values of y. Each intersection between a line and 
the curve corresponds to one image. If y > y r (bottom line), 
the source is outside the caustic circle, and only one image 
appears. If y = y r the source is on the caustic circle and 
a second image appears on the radial critical circle, at x = 
— x T - For y < y r , the source is inside the caustic, and the 
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Figure 4. Images of a circular source. Each pair of panels shows the source plane in the left panel, with the caustic, and the image 
plane in the right panel, with the radial (inner) and tangential (outer) critical circles. The position y of the source on the source plane 
is indicated. We used k c = 4.014, and a source of diameter Ay = 1. 



second image splits into two images. Finally, for y = (top 
curve), the central image is located at x = 0, and the two 
outer images are located on the tangential circle, at x — ±xt. 
Actually, these two images merge to form an Einstein ring. 
The slope of the curve a(x) versus x is equal to k c at x = 0. 
It is clear from Figure 5 that if we lower k c below 1, any 
x — y versus x line will intersect the curve only once; multiple 
images cannot occur if k c < 1. In Figure 6, we plotted the 
solution of the lens equation for various values of k c . The 
bifurcation from 1 to 3 images is clearly visible on all panels 
with k c > 1. 

In Figure 7, we plot the separation between the two 
outer images as a function of the source location, for various 
values of k c . The plot only extends to y/y r = 1, since larger 



values of y only produce one image. The separation is fairly 
insensitive to the source location, and stays within ~ 15% 
of the Einstein ring diameter Ax — 2xt for all values of 
K c considered. This is particularly convenient for theoretical 
studies, when the actual source location can be ignored. 

4.2 Magnification and brightness ratios 

The magnification of an image located at position x on the 
lens plane is given by 

d / m \ 1 — 1 

1 



(' 7) S ™ 
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Figure 5. Plot of ct(x) (solid curve) and x — y (dotted lines) 
versus x, for a TIS with ft c = 5.005 and 4 particular values of y: 
y = (top line), y = 4, y = y r = 7.515, and y = 12 (bottom line). 
Images are located at values of x corresponding to intersections 
between the lines and the curve. Particular solutions, correspond- 
ing to images located on critical curves, are indicated by arrows. 
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Figure 6. Location x of the image(s) on the image plane versus 
source location y on the source plane, for various values of n c . 




Figure 7. Separation Ax between the two outer images, in units 
of the tangential critical radius xt, versus source location y in 
units of the caustic radius y r . The solid curves correspond to 
various values of indicated. The dotted line at Ax = 2xt 

indicates the diameter of the Einstein ring. 
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Figure 8. Top four panels: Total magnification ^ttot versus source 
position y, for four different values of k c . Bottom panels: Bright- 
ness ratios /ii/fJ.3 and 112/^3 versus source position in units of 
caustic radius, y/y r - The number next to each curve gives the 
value of K c . 



(SEF, eq. [8.17]). We computed the total magnification and 
brightness ratios of images, for the four particular cases 
k c = 1.2, 2.5, 5.0, and 10.0. The results are shown in Fig- 
ure 8. The top four panels show the total magnification as a 
function of the source location y. As y decreases, the magni- 
fication slowly increases, until the sources reaches the radial 
caustic y = y r - At that moment, a second image, with infi- 



© 2002 RAS, MNRAS 000, [|-|l| 



8 Hugo Martel and Paul R. Shapiro 



nite magnification (1 — d(m/x)/dx = in eq. |24|) appears 
on the radial critical curve. As y keeps decreasing, that sec- 
ond image splits into two images, and the total magnification 
becomes finite again, until the source reaches y = 0, and an 
Einstein ring with infinite magnification (1 — m/x 2 — in 
eq. appears on the tangential critical curve. Of course, 
these infinite magnifications are not physical, since they can 
only occur for point sources. 

The total magnification is always larger than unity, and 
always larger when 3 images are present. Figure 8 shows that 
the total magnification decreases with increasing k c , which 
seems counter-intuitive. Notice, however, that the top four 
panels in Figure 8 are plotted for different ranges of y. At a 
fixed y, the total magnification always increases with k c . 

The bottom panels of Figure 8 shows the brightness 
ratios /ii//i3, and ^2/^3 1 where Hi is the magnification of 
image i, as a function of source position. By convention, 
image 1 is the one between the tangential and radial critical 
curves, image 2 is the one inside the radial critical curve, 
and image 3 is the one outside the tangential critical curve 
(see Fig. 4; from left to right, the three images are image 1, 
2, and 3). For most values of k c , image 3 is the brightest, 
image 1 is fainter, but comparable, and image 2 (the central 
one) is much fainter. This is a generic properties of most 
axially symmetric lenses with a central core. However, we 
find a different behavior at small values of k c . Not only can 
image 1 become brighter than image 3, but for a source 
located at y < y r , just inside the caustic, even image 2 can 
become brighter than image 3. 



4.3 Shear 

The total shear 7(2:) of an image located at position x is 
given by 



7 = 



m(x) 



k(x) 



(25) 



(SEF, eq. [8.15]), where k(x) = a(x)/ 'cr cr it- We substitute 
equations (M and (|l7|) into equation (|25j), and get 



7 = 



abKc 



Ab-Ba 



a + x 



A 



2\l/2 1 2B r. 2 



B 



(o 2 +a; 2 )V2 (&2 + 3.2)1/2 



(26) 



[This result was also obtained by Natarajan & Lynden-Bell 



(1997).] Figure 9 shows 7 versus source position y for various 
values of Kc- Notice that the indices 1, 2, 3 refer to the 
individual images, and not the the shear components. The 
shear increase with k c , as expected. At large values of « c , 
71 increase and 73 decrease as y increase, while at small k c 
this trend is reversed. 72 always increase with y. In the limit 
y — > 0, we get 71 = 73 and 72 = 0; Images 1 and 3 become 
identical, while image 2 becomes circular. 



4.4 Time delay 

For axially symmetric lenses, the deflection potential tp(x) 
is defined by a = dip/dx. From equation (|l4|), we get 

2abn c 



ij){x) = 



Ab-Ba 



[A{a 



2 . 2\l/2 
+ X ) 1 



Aa In [a + ( 



2 . 2x1/21 

a +x ) ' 
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Figure 9. Total shear 7 versus source position y, for four different 
values of k c . The indices identify the various images. 



-B(b 2 + x 2 ) 1/2 + Bb\n [b + (6 2 + x 2 ) 1/2 ] j 



(27) 



The time delay At between two images located at x% and 
Xj, of a source located at y, is given by 

rlD s 



At(y) 



■(1 + zl) [(f>(xi,y) - <t>(x},y)] 



(28) 



cD L D L s 

(SEF, eq. [5.44]), where zl is the redshift of the lens, and 
<j>(x, y) is the Fermat potential, defined by 



fay) = \i x - vf - i>( x ) 



(29) 



In equation (pq), the quantities Ds, Dl, Dls, and zt de- 
pend on the redshifts of the source and lens, and the back- 
ground cosmological model. Only the quantity in square 
bracket depends on the intrinsic properties of our lens model. 
We define 



(Xi,y) - <t>{xj,y) . 



(30) 



Figure 10 shows the time delays as functions of the source 
location for the cases k c = 1.2, 2.5, 5.0, and 10.0. In all 
cases, we have T23 > T13 ^ 0. Hence, the light from the 
three images always reach the observer in the same order, 
first image 3, then image 1, and finally image 2. As y — » 0, 
T13 — > because images 1 and 3 move toward the tangential 
critical circle. Both T13 and T23 increase with y, except for 
small values of k c , for which T23 reaches a maximum at some 
value y < y r , and then decreases slightly. At y = y r , we have 
T13 = T23 because images 1 and 2 merge on the radial critical 
curve. 

An interesting quantity is r max , the maximum time de- 
lay that a lens can produce. For simplicity, we set r max = 
r(y = y r ), which, according to Figure 10, is correct for large 
k c , and a fairly good approximation at small k c . In Fig- 
ure 11, we plot T max as a function of k c . The solution has 
the following asymptotic behavior, 

If: 



8.3(k c 



K c » 1. 



At small k c (but still larger than unity), r n 
value given by equation (fjlj). 



(31) 
exceeds the 
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Figure 10. Time delays T12 and T13 versus source position y in 
units of the caustic radius y r , for four different values of k c . 




log (/c c -l) 



Figure 11. Maximum time delays r max versus re c , on linear and 
logarithmic scales (solid curves). The dotted curves show the fit- 
ting formula given by equation fell). 



5 DISCUSSION 

So far, we have focused on the intrinsic, scale-free proper- 
ties of the lenses, without considering the dependence on 
the cosmological model. All results have been expressed in 
terms of two dimensionless parameters, k c and y. The cos- 
mological models enter the picture when one tries to deter- 
mine the typical values and distributions of these parame- 
ters. The value of k c depends on the critical surface density 
o"crit (eq. Jll]]), which is a function of the lens redshift, source 
redshift, and cosmological parameters. In addition, the TIS 



profile parameters po and ro are functions of the mass of the 
object and its epoch of collapse, which also depend upon the 
cosmological parameters. Finally, the distribution of source 
locations y will be related to the number densities of sources 
and lenses, which themselves depend upon the cosmological 
model. 

In this section, we will estimate the typical values of k c 
for TIS halos in a CDM universe. We consider the currently- 
favored ACDM model with density parameter Qo = 0.3, 
cosmological constant A = 0.7, and Hubble constant Ho = 
70kms _1 Mpc~* (or h = 0.7), and we shall assume for il- 
lustrative purposes that the source is located at redshift 
zs = 3. The critical surface density given by equation (jll]) 
depends upon the lens redshift zl, and has a minimum value 
at zl = 0.612, given by 
-2 



(a c 



0.386 gem" 



(32) 



The surface density of the TIS is computed by plugging the 
expressions for po and ro provided in Paper II into equa- 
tion (^). The central value of the surface density is given 
by 



cr(£ = 0) = irp r 



(a b) 
= 0.0400[F(z coll )] 2 ft 4/3 



M 



10 12 M o 



1/3 



gem 



(33) 



In equation (|33J), the function F(z co ll) is a scaling factor that 
is used to express the TIS solution for arbitrary cosmologi- 
cal models in terms of the solution for an Einstein-de Sitter 
model. We will further assume that the lens located at red- 
shift zl collapsed at redshift 2 co ii = zl- 

Let us consider the illustrative case of a lens at zl = 
0.612. Using the expressions from Paper II, we get F(z co n) = 
1.1644. Equation (p3J), with h — 0.7, reduces to 

1/3 

gem" 2 . 



<t(£ = 0) = 0.0337 



Mo 



\ 10 12 M 



Taking the ratio of equations 

1/3 



and (p2|), we get 



k c = 0.0874 



M 



10 12 M Q 



(34) 



(35) 



Since the phenomena of strong lensing (e.g. multiple images, 
arcs, . . .) requires k c > 1, equation (J33) indicates that, for 
our illustrative choice of (zl,zs) ~ (0.612,3) in a ACDM 
universe, strong lensing requires lensing halos as massive as 
M > 1.5 x 10 15 M Q . We can express this more generally 
by evaluating equations ([ll]) and (|33j) at different values 
of zl = Zcoii for a given source redshift zs, as shown in 
Figure 12. The quantity k c (M/10 15 M q )- 1/3 has a maximum 
value given by 



k c M, 



-1/3 



(z s = 3) = 1.1686, 



(36) 



for the ACDM universe (which occurs at zl ~ 1.320), where 
M15 = M/10 15 M Q . The condition k c > 1, therefore, re- 
quires 



M15 > 0.6266 . 



(37) 



To see how likely this is, let us focus on our illustrative 
example of (zl, zs) = (0.612, 3). For the cosmological model 
we consider, a 1-er density fluctuation collapsing at redshift 
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10 12 



Figure 12. The central convergence n c [divided by the factor 
(M /1O 15 M0) 1//3 which contains the dependence on halo mass 
M] of TIS halos in the ACDM universe versus lens redshift z^, 
for the illustrative case in which source is located at redshift 
25 = 3. The horizontal dotted line indicates the minimum value 
(k c = 1) required for strong lensing by cluster-mass halos of 
1O 15 M0. (Such halos can only produce strong lensing for re c > 1, 
i.e. 0.777 < z L < 1.914.) 



Zcoll = 0.612 has a mass of about 1.5 x 1O 12 M0 (the precise 
value depends on the details of the power spectrum and its 
normalization). Such "typical" objects will not be capable 
of producing multiple images of a source at redshift zs = 3, 
since the resulting value n c — 0.100 is smaller than unity. 
This simply indicates that multiple images are not produced 
by typical objects, which is certainly consistent with the fact 
that fewer than 100 multiple-image systems have been ob- 
served. Increasing k c above unity would require an object 
of mass M ~ 1.5 x 1O 15 M0, a thousand times the mass of 
a typical object at that redshift. Clusters of this mass are 
rare but do exist. We can make a simple estimate of how 
atypical such a massive object is. Over most of the mass 
range of cosmological interest (from small galaxies to clus- 
ters of galaxies) the CDM power spectrum can be roughly 
approximated by a power law P(k) oc k", where k is the 
wavenumber and n m —2. The rms density fluctuation <5 rms 
is then given by <5 rms « k 3/2 P 1/2 (k) oc k 1/2 . At a given 
redshift, different values of the wavenumber k correspond 
to different mass scales M according to Af oc k~ 3 . The rela- 
tion between rms density fluctuation and mass scale at fixed 
epoch is therefore approximated by 



oc Af 



-1/6 



(38) 



Increasing the mass by a factor of 1000 therefore reduces 
<5rms by a factor of 1000 1 / 6 w 3. Because of the reduction in 
<5rms, a 1-er fluctuation (<5 = <5 rms ) at this higher mass will no 
longer collapse by this redshift zl = 0.612, but a 3-a fluc- 
tuation (8 = 3<5 rms ) will. Such fluctuations are rare, but not 
vanishingly rare. In Gaussian statistics, the probability that 
a randomly located point in space is inside a 3-cr density 
fluctuation (i.e. 5 ^ 35 rms ) is 1/384. Of course, whether any 
halo produced by such a fluctuation will actually produce 



Figure 13. Ratio of the radius r i ng of an Einstein ring produced 
by a TIS halo to the radius $e of an Einstein ring produced by a 
Schwarzschild lens of the same mass (the Einstein radius), versus 



multiple images will depend on the location of the sources. 
We carried this simple calculation for illustration purpose 
only. In a future paper, we will present a detailed calcula- 
tion of the expected frequency of multiple image systems for 
comparison with the statistics of observed lensing. 

In the TIS model, there is a relationship between the 
parameters po and ro, and the total mass M of the halo, 
given by 



M = 4-irporo Aft 



(39) 



where M t = 61.485 (Paper I). Using this relation, we can 
directly estimate the separations between multiple images. 
The Einstein radius of a lens of mass M is defined as 
the angular radius 9e of an Einstein ring produced by a 
Schwarzschild lensfj] of the same mass, 

e B = (^^A y2 (40) 



D L D S 



(SEF, eq. [2.6a]). This Einstein radius is often used to es- 
timate the characteristic scale of image features caused by 
strong lensing (e.g. ring radius, radial location of arcs, im- 
age separations) and to estimate the size of the region within 
which the mass which is responsible for that strong lensing 
must be concentrated. Since lensing halos are not actually 
point masses, however, the angular radius ^ing of the actual 
Einstein ring which results if the source is located along the 
line of sight through the lens center will usually differ from 
the Einstein radius 6e, assuming that the lens mass distri- 
bution is actually capable of producing a ring. The angular 
radius # r m g of the actual Einstein ring produced by a TIS is 
given by 



_ r x t 



(41) 



where Xt is the solution of equation (|l8|). We take the ratio of 



A point mass. 
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equations ( [ioj ) and (^), and use equation (|ll]) to eliminate 
the angular diameter distances. We get, after some algebra, 



■Ting _ /7T\ J 

0E Vo 6 



M t - 1/2 K~ 1/2 x t 



= 0.1565«: c 1 ^ 2 x t . 



(42) 

In Figure 13, we plot 8 T i ng /0E as a function of k c . A ring 
can be produced only if k c > 1. The ratio of the actual 
ring radius to the Einstein radius monotonically increases 
with k c , and exceeds unity for n c > 9.096. Hence, the ring 
produced by a TIS can be either larger or smaller than the 
ring produced by a Schwarzschild lens of the same mass, de- 
pending on the value of k c . In the limit of large k c , we have 
Otiug/OE ~ 0.2563k;;!/ 2 according to equations ( pi] ) and (53). 
We saw in §4.1 that in cases of multiple images, the sepa- 
ration Ax between the outermost images is always of order 
2xt- Hence, 8 T i ng can be reinterpreted as being one half the 
separation between images. 



6 SUMMARY AND CONCLUSION 

We have derived the lensing properties of cosmological ha- 
los described by the Truncated Isothermal Sphere model. 
The solutions depend on the background cosmological model 
through the critical surface density (7 C rit, which is a function 
of the cosmological parameters and the source and lens red- 
shifts, and the TIS parameters po and ro, which are functions 
of the mass and collapse redshift of the halo, and the cos- 
mological parameters. By expressing the surface density of 
the halo in units of <r C rit and the distances in units of ro, all 
explicit dependences on the cosmological model disappear, 
and the solutions are entirely expressible in terms of two di- 
mensionless parameters, the central convergence k c and the 
scaled position y of the source. We have computed solutions, 
and we provide either analytical expressions or numerical 
fits, for the critical curves and caustics, the image separa- 
tions, the magnification and brightness ratios, the shear, and 
the time delay. Lensing of a point source by a TIS produces 
either one or three images, depending on whether the source 
is located outside or inside the radial caustic. When three 
images are produced, the central one is usually very faint, 
being highly demagnified. Degenerate image configurations 
occur when an extended source overlaps a caustic. Two im- 
ages are produced when the source overlaps the radial caus- 
tic, while an Einstein ring with a central spot is produced 
when the source overlaps the tangential caustic, which is 
a single point located at y t = 0. These degenerate cases 
correspond to maxima of the total magnification, which di- 
verges as the source size goes to zero. When three images 
are pro duced, the angular separation between the two out- 



ermost images depends strongly on k c , but only weakly on 
the source location. 

The lens properties are often qualitatively different at 
small and large n c . For instance, at small k c , k c > 1, the 
image separation Aa; decreases as source position y increases 
(i.e. as the projected separation between the source position 
and the lens center increases) (Fig. 7), the brightness ratio 
Mi / M3 increases with y (Fig. 8) , the shear 73 of the outermost 
image decreases with y, and the time delay T23 decreases 
with y for y < y r (Fig. 10), while these quantities behave 



imaging can only occur when the central surface density cr(0) 
exceeds the critical density a CI i t (or equivalently k c > 1). 
Since a(x) is a decreasing function of x, there is a natural 
scale in the system: the position :r cr it on the lens plane where 
c(:Tcrit) = o- cr it. If the lens profile was scale-free, as in the 
cases of a Schwarzschild lens or a singular isothermal sphere, 
Tcrit would be the only length scale in the problem, and the 
properties of the lens would be self-similar. But the TIS 
has a characteristic length scale, the core radius ro, and the 
existence of this second length scale prevents the solutions 
from being self-similar. 

This paper focused on the intrinsic properties of indi- 
vidual lenses described by the TIS model. It provides all 
the necessary formula one needs to study gravitational lens- 
ing by TIS halos in specific cosmological models. We will 
present such studies in forthcoming papers. As an illustra- 
tion here, we applied the TIS model to the currently-favored 
ACDM universe, to calculate the central convergence k c ex- 
pected for TIS halos of different masses and collapse epochs. 
We found that high-redshift sources (e.g. zs ~ 3) will be 
strongly lensed by TIS halos (i.e. k c > 1) only for cluster- 
mass halos. We also calculated the characteristic angular 
scale of image features produced by strong lensing by TIS 
halos relative to the Einstein radius 8e of a lens with the 
same total mass, for comparison with the results for other 
lensing halo mass profiles and with observed lensing systems. 
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